Magnetic interactions in EuTe epitaxial layers and EuTe/PbTe superlattices 
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The magnetic properties of antiferromagnetic (AFM) EuTe epitaxial layers and short period 
EuTe/PbTe superlattices (SLs), grown by molecular beam epitaxy on (111) BaF2 substrates, were 
studied by magnetization and neutron diffraction measurements. Considerable changes of the Neel 
temperature as a function of the EuTe layer thickness as well as of the strain state were found. A 
mean field model, taking into account the variation of the exchange constants with the strain-induced 
lattice distortions, and the nearest neighbor environment of a Eu atoms, was developed to explain 
the observed Tn changes in wide range of samples. Pronounced interlayer magnetic correlations have 
been revealed by neutron diffraction in EuTe/PbTe SLs with PbTe spacer thickness up to 60 A. 
The observed diffraction spectra were analyzed, in a kinematical approximation, assuming partial 
interlayer correlations characterized by an appropriate correlation parameter. The formation of 
interlayer correlations between the AFM EuTe layers across the nonmagnetic PbTe spacer was 
explained within a framework of a tight-binding model. In this model, the interlayer coupling stems 
from the dependence of the total electronic energy of the EuTe/PbTe SL on the spin configurations 
in adjacent EuTe layers. The influence of the EuTe and PbTe layer thickness fluctuations, inherent 
in the epitaxial growth process, on magnetic properties and interlayer coupling is discussed. 

PACS numbers: 75.70.Ak, 75.25. -hz, 68.65.Cd 



I. INTRODUCTION 

In the last decade magnetic multilayer systems and the 
giant magnetoresistance resulting from interlayer cou- 
plings have been receiving considerable interest in both 
applied and fundamental scientific research. Interlayer 
exchange couplings in multilayers and superlattices (SLs) 
have been observed in a large variety of structures com- 
posed of metallic ferromagnetic (FM) layers alternating 
with nonmagnetic metallio>i«^, as well as nonmetallioi^ 
spacer layers. These observations have stimulated exten- 
sive theoretical studies that have resulted in a number of 
different models for the mechanism of interlayer coupling 
such as the RKKY model, the free-electron model, and 
several others. The most complete theory unifying all 
previous approaches has been devised by Bruno^. How- 
ever, neither interlayer coupling in systems composed of 
two nonmetalUc materials, nor mechanisms that might 
give rise to coupling between antiferromagnetic (AFM) 
layers have been considered in these works. 

Yet, neutron diffraction data for three different SL 
systems composed of AFM and nonmagnetic semicon- 
ducting materials, reported in the mid-ninetie o^'7i^ii?i^ft , 
revealed the existence of pronounced interlayer corre- 
lations between the AFM blocks. Also recently, cou- 
pling between FM-semiconductor layers has been found 
in EuS/PbS SLaii. In these all-semiconductor systems, 
the carrier concentration is far too low to support any sig- 



nificant RKKY interactions; in addition, the AFM layers 
do not carry a net magnetic moment. Thus, the two main 
ingredients which were believed to play a crucial role in 
interlayer coupling - mobile carriers and layer magnetiza- 
tion - are absent in these cases. These results have clearly 
demonstrated that the magnetic interlayer coupling is not 
restricted to structures containing FM metallic compo- 
nents. The proper understanding of correlations between 
the AFM semiconductor layers may shed new Hght on the 
interlayer coupling mechanisms. 



In this paper we present our systematic, experimen- 
tal and theoretical, studies of EuTe epitaxial films and 
short period [(EuTe)m|(PbTe)„]Ar superlattices. In Sec- 
tion^we describe the basic properties of the constituent 
materials, sample preparation and the experimental tech- 
niques used. The effects of the finite thickness and the 
strain on magnetic properties of individual layers are 
treated in Section UTTl Section ITvl is devoted to the inter- 
layer coupling determined by neutron diffraction mea- 
surements. In its first three subsections, the neutron 
data for a series of samples studied in zero, intermedi- 
ate and high magnetic in-plain field, as well as the effect 
of cooling in external low magnetic field are presented. 
In the last subsection we discuss the results of a the- 
oretical model for the interlayer coupling in a perfect 
AFM/nonmagnetic semiconductor SL and we compare 
the experimental results with the model predictions. 
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FIG. 1: Chemical and magnetic unit cell of EuTe with its 
type II AFM structure consisting of ferromagnetically ordered 
(111) planes and antiferromagnetic spin sequences along the 
[111] directions. 



II. CONSTITUENT MATERIALS, SAMPLE 
PREPARATION AND EXPERIMENTAL 
TECHNIQUES 

Bulk EuTe is a classical Heisenberg antiferromagnet 
with a Neel temperature Tn = 9.6 Kli. It crystallizes in 
the NaCl structure with a = 6.598 A. The Eu^+ ions 
with S = 7/2 and L = form an FCC spin lattice 
with dominant AFM next-nearest neighbor interactions, 
and weaker FM interaction between the nearest neigh- 
bors ( J2/fcB = -0.23 K, Ji/fce = 0.11 K, respectivelj*i^) . 
Such a Ji combination leads to the Type II AFM order- 
ing below Tn, in which the spins are ferromagneticahy 
ordered in the (111) lattice planes, and the neighboring 
planes are coupled antiferromagnetically to one anothesi^ 
(see Fig.^. EuTe is a wide gap 2.5 eV) semiconductor 
with the 4/ levels situated about 2 eV below the conduc- 
tion band edge^^. The diamagnetic constituent, PbTe, 
is a narrow gap (~ 0.19 eV) semiconductor, which also 
crystallizes in the NaCl structure and has a bulk lattice 
constant of 6.462 A. This yields a close lattice-matching 
to EuTe, with a mismatch of only 2.1 % in the lattice 
constants. 

The EuTe/PbTe SL samples were grown by molec- 
ular beam epitaxy on (111) oriented BaF2 substrates 
as described in detail in Ref. Different 
[(EuTe)m|(PbTe)„]Ar SL stacks with m (varying from 2 to 
10) monolayers of EuTe alternating with n (from 5 to 30) 
monolayers of PbTe, were deposited on 1 - 3 fim PbTe 
buffer layers. In the paper the EuTe/PbTe SL with m 
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FIG. 2: (a) Evolution of the (^^|) magnetic Bragg peak from 
3 /im EuTe epilayer with the temperature. The asymmetry of 
the peak indicates that the relaxation of the strain in the EuTe 
film is not complete. (6) Integrated intensity of the (HI) and 

(HI) peaks vs. temperature for the same epilayer. (c) The 
magnetic moment vs. temperature of 3 /im EuTe epilayer 
measured by SQUID in applied magnetic field of 1 Tesla. 



monolayers of EuTe and n monolayers of PbTe in the SL 
period (bilayer) is often denoted by an abbreviated sym- 
bol {m/n). The thickness of one monolayer is 3.81 A for 
EuTe and 3.73 A for PbTe. To obtain sufficient neutron 
scattering intensity, the number of periods N was several 
hundred in all cases. The layer thicknesses determined by 
high-resolution x-ray diffraction agreed with the nominal 
thicknesses within ±0.5 monolayer. The electron concen- 
tration in the PbTe layers was ~ 10^ ''cm"'^, i.e., many 
orders of magnitude lower than in metals, and the EuTe 
layers were semi-insulating. 

The neutron experiments were performed at the NIST 
Center for Neutron Research. The instruments used 
were BT-2 and BT-9 triple-axis spectrometers set to 
elastic diffraction mode, with a pyrolitic graphite (PG) 
monochromator and analyzer, and a 5 cm PG filter in the 
incident beam. The wavelength used was A = 2.35 A and 
the angular collimation was 40 minutes of arc through- 
out. Additionally, a number of diffraction experiments 
were carried out on the NG-1 refiectometer operated at 
neutron wavelength A = 4.75 A. The latter instrument 
yielded a high intensity, high resolution spectra with a 
neghgible instrumental broadening of the SL diffraction 
lines. All the magnetic diffraction spectra reported here 
have been measured at 4.3 K. 

The dominant feature in diffraction spectra from Type 
II antiferromagnets is a strong refiection at the (555) 



3 



position. Pronounced maxima centered at the (555) 
position were observed in all EuTe and EuTe/PbTe SL 
samples cooled below the Neel temperature, including 
those in which the EuTe layer thickness was as small as 
2 monolayers. The magnetic origin of the AFM SL peaks 
was confirmed in several ways. First, the SL reflections 
also appear at other Q-space points characteristic for the 
AFM II structure, their intensities being consistent with 
the Eu2+ magnetic formfactor. Secondly, the scattered 
intensity shows the typical temperature behavior, closely 
following the squared Brillouin function for a S ~ 7/2 
system below Tn. Finally, we have performed a number 
of diffraction measurements in external magnetic fields 
that convincingly prove the maxima, we have investi- 
gated, originate from the ordering of the Eu magnetic 
moments and not from other effects that may potentially 
produce peaks at the same reciprocal lattice points (such 
as, e.g., chemical ordering). 

The maxima observed by us show that the Type II 
AFM ordering was preserved in all EuTe layers and 
EuTe/PbTe SLs studied. Shown in Fig. Illa,^) is the 
temperature dependence of the magnetic peaks intensity 
for 3 /xm EuTe epilayer sample. In Fig.[2lc) the temper- 
ature dependence of the magnetic moment of another 3 
/im EuTe epilayer measured by a standard SQUID tech- 
nique is presentedi^. The critical temperatures deter- 
mined by neutron diffraction and magnetization mea- 
surements agree within the experimental errors. The ob- 
tained value of Tn — 10.4 ±0.05 K is slightly higher than 
9.6 K found in bulk EuTe, due to the strain introduced to 
the epitaxial film by the BaF2 substrate and PbTe buffer 
layer. The distorted, non-gaussian profile of the (555) 
magnetic Bragg peak in Fig. [2ji clearly points out to the 
nonuniform lattice distortions present in the sample, the 
closer the portions of the EuTe film to the substrate the 
stronger the deformations of the EuTe lattice. Thus in a 
sense, the 3 /xm layer constitutes only a semi-bulk sam- 
ple. The influence of strain on magnetic properties of 
EuTe layers will be discussed in detail in the following 
section. 



III. STRAIN AND FINITE SIZE EFFECTS 

A. T-domain structure 

In a perfect FCC lattice there are four symmetry- 
equivalent Type II AFM arrangements in which the FM 
spin sheets form on the (111), (Til), (iTl), or (111) 
plane families (see Fig. El . These four configurations are 
usually referred to as "T-domains". In a macroscopic, 
strain-free EuTe crystal cooled through the Neel point at 
Hext = all four T-domain states would be populated, 
giving rise to magnetic peaks at (^55), (555)) (555) 
and (5^5) reflection points. However, in the EuTe epi- 
layer and SL samples only the (555) maximum is ob- 
served in the neutron diffraction spectra; no detectable 
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FIG. 3: The four symmetry equivalent Type II AFM arrange- 
ments of the {111} ferromagnetic spin sheets in EuTe. 



TABLE L Neighborhood spin configuration in an EuTe layer 
for different T-domain types in which the ferromagnetic spin 
sheets are either parallel ('in-plane') or inclined ('oblique') to 
the (111) epitaxial growth plane. 



Monolayer 
number: 


Domain type: 


'in-plane' 


'oblique' 


NNs 


NNNs 


NNs 


NNNs 


i - 1 


3(i) 


3(i) 


2(T)+ l(i) 


3(i) 


i 


6(T) 


none 


2(T)+ 4(i) 


none 


i + 1 


3(i) 


3(i) 


2(T)+ l(i) 


3(i) 



magnetic scattering was ever found at the other three re- 
flection sites in any of the investigated specimens. This 
means that in the layer systems only a single T-domain 
state forms - the one in which the FM spin sheets are 
parallel to the EuTe layers - whereas the other three 
'oblique' configurations never occur. 

The observed preference in the T-domain formation 
can be explained by simple energy arguments. In the 
Type II AFM structure a spin residing in a given FM 
sheet has six FM nearest neighbors (NNs) located in the 
same sheet, and three AFM ('frustrated') NNs in each ad- 
jacent sheet. The AFM next-nearest neighbors (NNNs) 
are also located in the adjacent sheets, three in each. 
Suppose that an (111) EuTe layer consists of m mono- 
layers. Consider, e.g., an 'up' (t) spin located in the i^^ 
monolayer. TableQlshows its neighbor environment (i.e., 
how many of its NNs and NNNs with parallel (t) and an- 
tiparallel (J,) spins are located in the (i — 1)*'^, (i)*'\ and 
the {i + monolayer) for the 'in-plane' and 'oblique' 
domain configurations. Because all NNNs in the Type 
II AFM structure are antiparallel, the total energy of 
the NNN interactions is the same in the 'in-plane' and 
the 'obHque' domains. The spins residing in monolayers 
with numbers 2 < i < m — 1 have a full set of six parallel 
and six antiparallel NNs, so that the total energy of the 
NN interactions for these spins is zero, regardless of the 
domain arrangement. However, the spins in the inter- 
face monolayers (i.e., those with i = 1 and i — m) have 
only nine magnetic NNs: 6(t)+3(i) for the 'in-plane' and 
4(T)+5(i) for the 'oblique' domains. Hence, due to the fi- 
nite thickness of the layer, the exchange energy of the NN 
spins is not equal to zero any more and different for the 
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'in-plane' and 'oblique' configurations, with correspond- 
ing average magnetic energy per Eu atom of : 



and 



^obl 



2C 



= 2C 



-3Ji 3(to-1] 



J2 



Ji 3(m - 1) 



J2 



(1) 



(2) 



where C = S{S + l)/3. Because Ji > 0, the 'in-plane' 
spin arrangement is the one that minimizes the total mag- 
netic energy. 

Another factor that favors the 'in-plane' domain con- 
figuration is strain. Since apbXe < lEuTe, the EuTe lattice 
is 'compressed' in the layer plane. Therefore the distance 
between the NN Eu ions located in the same (111) 
monolayer shortens, while, due to the lattice reaction to 
the strain the distance between the NNs residing 

in the adjacent (111) monolayers increases. The dj^j^ 
and values may differ from the bulk NN distance 
d^^ by as much as ±2%. As shown by Goncharenko 
and Mirabeauii from neutron measurements under high 
hydrostatic pressure, the Ji exchange constants in Eu 
chalcogenides generally increases with the decrease of the 
ion-ion separation. Hence, one can expect the jf for the 
'in-plane' NNs to be higher, and the for the 'out- 
of-plane' NNs lower than J}"""^. However, for EuTe the 
rate of Ji change with c^nn in the antiferromagnetic state 
could not be measured in Ref. because the contribu- 
tion of Ji to the Neel temperature cancels in the cubic 
AFM state, i.e., Ji is not accessible by experiment under 
hydrostatic pressure. With respect to the NN exchange, 
J2 appears to be only weaklydependent on the ion-ion 
separation, as shown in Ref. [13 . 

Taking into account both the finite thickness and strain 
factors, we obtain the average magnetic energy per Eu 
spin in EuTe layer consisting of m monolayers for the 
two different domain arrangements in the form: 



2C 



and 



,obl 



2C 



AJi 



m 



3(m- 1) 



J2 



(3) 



(4) 



where AJi = jj' — J^. 

The first right-side term, which does not depend on the 
layer thickness, refiects the effect of strain in the observed 
domain type preference. The second term represents the 
finite thickness effect as already given in Eq. Q and l|2l . 
In the thick layer limit (m ^ 00) the energy becomes 

fi" = 2C[-3A Ji + 3 J2] 
e°" = 2C[A Ji + 3 J2] 



which for unstrained samples (A Ji = 0) leads to the bulk 
value £^1""= = 6CJ2 for both arrangements. For small m 
values, however, the difference becomes significant (e.g., 
for m — 5 and m — i the energy per spin for the in- 
plane domain arrangement is lower, respectively, by 8% 
and 35% than for the obHque ones). 

For the in-plane compressed EuTe layers we introduce 
a parameter which describes the ratio of the elongation 
of the 'out-of-plane' NNs bonds to the shortening of the 
distance between the 'in-plane' NNs: k = Ad^^ / Ad^j^-^ . 
This parameter can be expressed in terms of the in-plane 
eii and out-of-plane e± strain, i.e.. 



i=y^(l + e||)2 + 2(1 + 6^)2 



LV3 



(5) 



For small strain values e, k is essentially constant and 
directly proportional to the e±/e\i ratio. With the value 
of the Poisson ratio for biaxially strained (111) EuTe lay- 
ers viii — 0.301 and the relation e±/e^\ — 2v/{v~ 1), we 
obtain /c « —0.24 for a compressive in-plane strain values 
of e|| < 2%. 

Using the k parameter and denoting by ^ the value of 
(SJi/Mnn), we can rewrite the magnetic energy given 
by Eq. (pj, in the form: 



-6C 



Ji (to - 1) 



J2 



{m{k + 1) - 1) 



erf] 



NN 



(6) 

In the EuTe/PbTe structures the compressive strain 
and finite thickness both appear to favor the 'in-plane' T- 
domain arrangement. It is interesting to note, however, 
that if the spacer material had a larger lattice constant 
than EuTe and produced a tensile strain in the EuTe 
layers, the strain would favor the 'obHque' arrangement 
and the two effects would compete with each other. The 
domain arrangement in such superlattices could be then 
tailored by manipulating the EuTe layer thickness. 



B. Changes of the Neel temperature 

The effect of the strain and the finite thickness in the 
EuTe/PbTe system is not only demonstrated by the pre- 
ferred in-plane spin alignment in the EuTe layers but 
also substantial shifts in the Neel temperature. This is 
clearly indicated by the measurements of the intensity 
of the magnetic diffraction signal as a function of 

temperature as shown in Fig. 31 for several different SL 
samples with different EuTe layer thicknesses to. The 
shapes of these curves were found to be in good agree- 
ment with the squared mean-field Brillouin magnetiza- 
tion function for S = 7/2 indicated by the soHd fines in 
Fig. 31 The transition temperature Tn was determined 
by fitting the function to the measured data, with Tn as 
an adjustable parameter. For most samples the experi- 
mental Tn values differs significantly from the bulk value 
of Tn = 9.6 K (see Fig.|l||. 
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FIG. 4: The intensity of the (t 



magnetic diffraction peak 



for different EuTe/PbTe superlattices. The solid lines corre- 
spond to the fit with the squared mean field Brillouin magne- 
tization function, with Tn as an adjustable parameter. 
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FIG. 5: Magnetic susceptibility of several EuTe/PbTe super- 
lattices measured with 10 Hz AC SQUID magnetometry. 



The SL samples were also investigated by magnetomet- 
ric methods^. Examples of the magnetic susceptibility 
vs. temperature dependence, which was measured using 
a 10 Hz AC SQUID magnetometer, are shown in Fig. El 
for samples with constant EuTe layer thickness of m = 5 
(top panel) and m = 10 (lower panel) but varying PbTe 
spacer thickness n. The observed x(^) characteristics 
for the SL specimens showed distinct maxima indicative 
for the AFM phase transition. The Neel temperatures 
obtained by this method were in good agreement (within 
±0.4 K) with the Tn values yielded by neutron diffraction 
experiments. 

The analysis of both, the neutron diffraction 



and magnetization data for a large number of 
[(EuTe)m|(PbTe)„]Ar specimens with different m and n 
values reveals certain distinct trends in the behavior. 
The experiments were carried out on series of SL samples 
for which: 

(i) either the strain in the EuTe layers was approxi- 
mately constant (by keeping the ratio of to to n 
fixed), see Figs. and Ela), or 

(ii) the EuTe layer thickness (to) was constant, and the 
PbTe spacer thickness (n) varied, see Figs. and 

mb). 

In the 'constant-strain' samples with very thin EuTe lay- 
ers (to = 2) the Tn was found to be considerably lower 
than the bulk value 9.6 K, but it increased with increas- 
ing TO. For TO « 5 it exceeded the T^""^, showing a ten- 
dency to level out at a significantly higher value of 12.6 
K (Fig. Ela)). For the structures with fixed EuTe layer 
thickness the Tn showed a clear growing tendency when 
the PbTe spacer thickness was increased. The larger 
value of n for constant to increases the strain in the mag- 
netic layers. The in-plane strain ey (or the in-plane lattice 
constant ay) within the EuTe layers was determined by 
x-ray diffraction. Plotted against ay, the Tn data from 
the sample series with constant to show approximately 
linear behavior (Fig.Et6)). 

The trends observed in the Tn behavior can be ex- 
plained on the grounds of the same simple model that has 
been used for explaining the preference in the T-domain 
formation. Adopting the mean-field theory approach, one 
can assume that the phase transition temperature is pro- 
portional to the effective field experienced by an 'average' 
spin at T = - in other words, to the average energy per 
spin e in the ground state, given by Eq. (EJ. Taking 
into consideration that Tn/Tj^""^ — e/e''""^, one obtains 
the expression for the Neel temperature of thick strained 
layers: 



N — Jn 



bulk 



1 



AJi 



(7) 



For small lattice distortions we assume that A Ji is pro- 
portional to the distortion parameter and hence, Tn 
should exhibit a linear dependence on the in-plane lat- 
tice constant. From the linear fit to the to = 10 data 
points in Fig. ^b) we obtain the rate of change of Tn 
with a||: dT^/da\\ = —33.8 K/A±5% and, correspond- 
ingly, ^ = 0.41 K/(AfcB). For thinner layers, when the 
—&Ji/m term may not be neglected, this linear charac- 
teristic should shift to the left. This is indeed consistent 
with the Tn vs. ay behavior observed for the m = 10 and 
m — 5 sample series. 

To obtain an expression for Tn as a function of the 
distortion parameter AdJI^jy and the layer thickness to, 
we use as before the Tn/Tj^""^ = e/e^uik relation with e 
given by Eq. JHJ: 



Tn 



T^bulk 



1- 



1 



1 



A 

\J2\ 



toIJqI 



(8) 
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FIG. 6: (a) The Neel temperature for the 'constant strain' 
samples (m:n=l:3) vs. EuTe layer thickness, (b) The strain 
dependence of the Neel temperature for the two families (5/n) 
and (10/n) of EuTe/PbTe superlattices. The solid lines are 
only a guide for the eye. 

From Eq. it follows that: 

(i) For samples with the same to, Tn exhibits a linear 
dependence on Adj^j^. 

(ii) The slope of the Tn vs. Ad|[j^ line is equal to 

^Ac?||jj^,(2to — 1)/to| J2I; it is weakly dependent on 
TO, except for very thin layers, where the change in 
the (2to — 1)/to factor becomes significant; 

(iii) With decreasing m, Tn shifts to lower values. The 
'zero-strain' value is given by Eq. Q by extrap- 
olation to Adjljj^ = 0). 

The qualitative predictions of the above simple model 
appear to be consistent with the Tn vs. Adjl^j^ mea- 
sured data for the m = 10 and to = 5 SLs series. The 
rate of Ji change yielded by the data appears to be ap- 
proximately 24% per 1% in the ion-ion distance change. 
This result supports the conclusions by Goncharenko and 
Mirabeauii concerning the Ji dependence on distance. 



The calculated 'zero-strain' Tn for m = 10 and to = 5 
is 9.1 K and 8.6 K, respectively, whereas the correspond- 
ing values extrapolated from the measured characteristics 
are 8.5 K and 7.0 K. From Eq. ijHJ one obtains that the 
lowest possible Neel temperature is 7.2 K (for unstrained 
layers with to = 2). Actually, experiments on samples 
with to = 2 yielded Tn values close to 5 K. 

This simple mean-field model correctly explains the 
qualitative behavior of Tn in layers with various thick- 
nesses and strain values. However, quantitatively the 
model appears to be less successful, especially for very 
thin layers, showing 10% to 20% discrepancy between 
the model and experimental Tn values. One possible 
reason of this discrepancy may be structural imperfec- 
tions that certainly exist in real superlattices. For in- 
stance, as indicated by the results of magnetization stud- 
ies of EuTe/PbTe(100) SLsiS and of a similar SL system 
EuS/PbSifii, even slight interdiffusion effects in the mag- 
netic/nonmagnetic interface regions may lead to observ- 
able decrease of the phase transition temperature due to 
the reduction of the number of NN and NNN exchange 
bonds between the Eu spins. 



C. S-domain structure and net magnetic moment 

In the model outlined above it was assumed that the 
Eu spins in each individual EuTe layer form a perfectly 
homogeneous Type II AFM order. The fact is, however, 
that in the (111) layer plane there are three (112) easy 
axes, 120° apart. This makes possible six microscopically 
inequivalent domain arrangements (usually referred as 
the S'-domainsSSi). It becomes a natural question whether 
each individual EuTe layer in the SL structure constitutes 
a single ^-domain, or does it consist of many smaller 
ones, in which the spins are oriented along different easy 
axes. Another important question concerns the mag- 
netic moment of the SL. The layers can be thought of 
as truly AFM only if to is an even number. For an odd 
TO, however, the layer as a whole, or the constituent S*- 
domains should posses an uncompensated moment - in 
other words, the layers become ferrimagnetic. Note, that 
the opposite spin configurations in successive EuTe layers 
could also lead to zero net magnetic moment of the entire 
SL, but this would require strong, perfect interlayer spin 
correlations. 

Information about the S'-domain structure and the net 
layer moment is important for understanding the inter- 
layer coupling effects seen in the EuTe/PbTe SLs. There 
is no direct method of visualizing domains buried in a 
SL structure or measuring their uncompensated moment. 
Yet, much insight into both these issues may be ob- 
tained from magnetization measurements and from neu- 
tron diffraction studies of spin rotation processes in a 
magnetic field applied parallel to the EuTe layers. 

In standard measurements using an unpolarized neu- 
tron beam the magnetic diffraction intensity is propor- 
tional to cos^ a, where a is the angle between the spins 
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FIG. 7: Directions of Eu spins Sj in zero external field (a), 
and in the applied field < 1 T (6). 



and the reflecting plane. Since in the EuTe/PbTe sys- 
tems the spins He in the (111) plane, for the {^W) re- 
flection a = 0. Any spin rotation in the (111) plane 
does not change this value, and thus the reflection in- 
tensity does not change. The information about spin 
orientation can be obtained from this reflection only by 
using polarized incident neutrons and polarization anal- 
ysis of the diffracted beam. However, as the intensity in 
polarized neutron experiments is typically about an or- 
der of magnitude or more lower than in measurements 
with unpolarized beam, such studies appear to be too 
time-consuming in the case of the EuTe/PbTe multilay- 
ers. Fortunately, the same information can be obtained 
by studying the (Hf) magnetic reflection, taking ad- 
vantage of the fact that the reflecting plane associated 
with it, (113), is nearly perpendicular to the (111) plane 
(arccos l/-\/33 — 80°), see Fig. In such experiments 
the external fleld i?ext is appHed parallel to the [110] 
axis (i.e., the axis of intersection of the (111) and (113) 

planes), and the (^^f) reflection intensity is measured 
vs. the fleld strength, as shown in Fig. |H1 

First, let us discuss the anticipated outcome of such an 
experiment based on the ideaHzed picture of a EuTe layer. 
Suppose that all possible S'-domain states are equally 
populated. Hence, when i?ext — 0, one-third of all spins 
are parallel to the [112] easy axis and make an angle 
a — 80° with the (113) plane. Two-third of the spins He 
along the [211] and the [121] easy axes (Fig.Ua)), and for 
them a = 29.5°. Thus, the observed reflection intensity 
is proportional to / cx ^ cos^ 80° | cos^ 29.5° = 0.515. 

When the fleld is turned on, the system reaction should 
depend very much on whether the number of spin mono- 
layers is even or odd. Consider flrst an even to, so that 
there is no uncompensated moment. In such a situation, 
the fleld tends to turn the spins toward an orientation 
perpendicular to /^ext- It does not affect the orientation 
of I of the spins along the [112] axis which are already 
perpendicular to i?ext- To turn the remaining 
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FIG. 8: Intensity of the (ii|) refiection vs. applied magnetic 
field for a bulk EuTe specimen and several EuTe/PbTe SLs. 



spins, the fleld has to overcome the anisotropy; the fleld 
strength needed for that is i?ext ~ VH^Hj, where Ifa 
is the anisotropy fleld, and Hj is the exchange fleld ex- 
pressed in magnetic units. Since in EuTe has been 
found to be about 12 G at 4.2 K, and Hj is about 3.5 T, 
the flop to the perpendicular position should occur be- 
fore i/ext reaches 650 G. After that, all spins make an 80° 
angle with the reflecting plane (see Fig.Effe)), and the re- 
flection intensity is / oc cos^ 80° = 0.03; in other words, 
it should drop to 5.9% of its zero-fleld value. With fur- 
ther iJext increase the intensity should not change much 
until reaching the Tesla region, where the Eu moments 
start inclining signiflcantly toward the fleld, leading to 
further suppression of AFM reflection. However, when 
the fleld is graduaHy decreased, the original zero-fleld 
intensity should not be restored, because aH the spins 
should remain 'locked' in the perpendicular position by 
the anisotropy fleld. 

This idealized model scenario changes quite dramati- 
caHy for an odd to. Now, the fleld tends to aHgn the un- 
compensated moment parallel to i?ext • The fleld strength 
needed to overcome the anisotropy is ~ mHa, i.e., < 100 
G for the to < 10 systems studied by us. In such i?ext, aH 
the spins get aligned parallel to the (113) plane, so now 
a = and / esc cosO = 1; in other words, in i?ext ~ 100 
G, there should be an increase of the (^^f ) reflection in- 
tensity by a factor of almost two. After that initial jump, 
the intensity should not change until the fleld reaches a 
value where inclining of the Eu moments from the mag- 
netic fleld lowers the system energy. After the fleld re- 
turns to zero, the spins should chose easy axes nearest 
to the fleld direction, so that the S'-domains correspond- 
ing to the [112] direction should not be repopulated. All 
spins now make a 29.5° angle with the reflecting plane, 
so that / oc cos 29.5° = 0.8 should be about 60% higher 
than originally. 

As illustrated in Fig. |H1 the signiflcant increase of the 
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reflection intensity, expected to occur in samples with 
odd TO, was never observed in any real EuTe/PbTe SL. 
A slight increase (~ 10%) was observed in some speci- 
mens - however, for samples with even m. Moreover, the 
observed /(77ext) characteristics are almost completely 
reversible - there is no indication for any hysteresis. Af- 
ter the field was increased and then returend back to 
zero, in all samples the original reflection intensity was 
restored within experimental error. 

All studied samples behaved essentially as in the 'even- 
TO scenario' (apart from reversibility), thus showing no 
uncompensated magnetic moment. However, the fleld 
needed to rotate all spins to the direction perpendicular 
to -ffext was always significantly higher than the expected 
value of i?ext = V Hg^Hj w 650 G, corresponding to the 
bulk value of Hg_ — 12 G. Moreover, the shape of the 
curves in Fig.|Slis not consistent with a single i/g value, 
but rather suggest that in each sample there is a sta- 
tistical distribution of the anisotropy fields. In fact, a 
satisfactory description of the observed curves was ob- 
tained by assuming a Gaussian distribution of the i?a 
values. The mean i/a values obtained from the fits for 
different samples varied from 50 to 200 G, but no system- 
atic trend has been found. To comment on the fact that 
no traces of any hysteresis were observed in the studied 
samples, we recall that the expected 'even-m scenario' 
irreversibility was deduced for idealized, perfect layers 
with the three-fold symmetry of anisotropy fields as in 
the bulk material. In real SLs, not only are the values of 
-ffa different, as obtained from the fits in Fig.|Hl but also 
distribution of their directions may deviate substantially 
from the <211> axes due to the infiuence of various types 
of defects and inhomogeneities. 

The absence of the ferrimagnetic properties in the 
SL samples with nominally odd number m of magnetic 
monolayers, seen in the (^^f) refiection intensity vs. 
field strength measurements, was further corroborated 
by the magnetization studies. In Fig.|2|the magnetic mo- 
ments for various samples, per SL period, per mm^, are 
plotted vs. the magnetic field up to 1 kG, which should 
align all the uncompensated spins parallel to the field. As 
shown in Fig.|2l for all SL samples, with both nominally 
even and odd to, the magnetic moment is several times 
smaller than the value 340 nanoemu one should observe 
from the 1 mm ^ of an uncompensated EuTe monolayer. 

The observed field dependence of the intensity of the 
(iM) psak, the absence of hysteresis and the low val- 
ues of the net magnetic moments at intermediate mag- 
netic fields all indicate that the vanishing of the magnetic 
moments in the SLs with nominally odd numbers of Eu 
monolayers can neither result from a random orientations 
of ferrimagnetic domains, nor from an antiparallel orien- 
tation of such domains in consecutive EuTe layers due 
to interlayer coupling. It must result from almost com- 
plete compensation of the magnetic moment within each 
domain. Such an unexpected compensation may be at- 
tributed to, e.g., a specific terrace structure of EuTe lay- 
ers with one monolayer steps as shown in Fig. EH The 
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FIG. 9: Magnetic moment per SL period, per mm^ for sev- 
eral EuTe/PbTe SLs with nominally odd and even number 
of magnetic monolayers in EuTe layer vs. applied magnetic 
field. 
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FIG. 10: Diagram of (5/4) EuTe/PbTe SL with conformally 
repeated terrace structure, that leads to a large reduction 
of the total magnetic moment of the individual odd-m EuTe 
layer as compared with an identical layer without such a step. 
For clarity, only the cations, Eu (solid circles) and Pb (open 
circles), are presented and the anion Te atoms are omitted. 



existence of one monolayer thick steps on the surface 
of MBE deposited EuTe layers has been confirmed by 
scanning tunneling microscopy studiesaa. The X -ray and 
neutron diffraction spectra, reported in the next section, 
prove the very high structural quality of our SLs, with- 
out traces of any significant interface roughness. They do 
not exclude, however, such terrace structures, the more 
so when the steps are conformally repeated over several 
SL periods. 
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FIG. 11: Possible spin configurations in EuTe/PbTe SLs. (a) 
Uncorrelated SL; directions of the monolayer magnetizations 
in consecutive EuTe layers change randomly giving rise to a 
single broad maximum in a neutron diffraction pattern, (b) 
Correlated SL; directions of the monolayer magnetizations in 
consecutive EuTe layers change in a regular way; in the case 
presented in the Figure, the orientations of the spins in all lay- 
ers are the same. Corresponding diffraction pattern exhibits 
a number of narrow fringes. 



IV. INTERLAYER COUPLING 

A. Neutron diffraction in zero magnetic field 

Magnetic neutron diflFraction is the only experimental 
tool capable of revealing the interlayer spin correlations 
in the case of AFM/nonmagnetic multilayers. The prin- 
ciple of the method is illustrated in Fig. Neutron 
diffraction scan along the [111] direction through the 
(^^^) reflection point from [(EuTe)io|(PbTe)3o]3oo speci- 
men with rather large thickness of the nonmagnetic PbTe 
spacers (rfpbTe = 112 A) is displayed in Fig. I12r a'). The 
large number of SL satellite peaks in the X-ray diffraction 
spectrum and the excellent agreement with dynamical 
simulations (Fis. \V2( b)) prove that the structural quality 
of the specimen is very high. In contrast to the multi- 
peaked X-ray pattern, the magnetic neutron diffraction 
spectrum of this SL has only the form of a single broad 
peak accompanied by two weak subsidiary side maxima. 
This profile shows a close similarity to the squared struc- 
ture factor of a single (EuTe)m layer, |FblP, known from 
the standard diffraction theory^J (also see the Appendix) . 
Such a spectrum shape produced by a multilayer struc- 
ture indicates the lack of coherence between the waves 
scattered by successive layers, meaning that the spin 
alignments in these layers are not correlated. However, 
when the PbTe spacer thickness decreases, the character 




en 

FIG. 12: (a) Neutron diffraction scans along the [111] direc- 
tion through the principal (^||) AFM reflection in a 10/30 
and 4/12 SL sample, showing pronounced satellite peaks for 
the latter one. The dotted line represents the expected shape 
of the spectrum for the fully correlated SL, with the instru- 
mental resolution taken into account. The dashed line shows 
the single layer magnetic structure factor. (6) High-resolution 
X-ray diffraction scan through the (222) reflection in the 
10/30 SL specimen (points: measured data, line: dynamical 
simulation). 



of the AFM reflections dramatically changes. As exem- 
plified in Fig. El for the 4/12 SL, a distinct pattern of 
narrower satellite peaks then emerges at regular intervals 
AQz equal to the spacing between the satellite peaks in 
the X-ray spectra. This clearly indicates the formation of 
magnetic interlayer correlations across the PbTe spacers. 
For dpbTe below ~ 60 A these magnetic satellites become 
the dominant part of the spectrum, as is also shown in 
Fig. El for a series of [(EuTe)5|(PbTe)„]3oo SL samples 
with varying PbTe spacer thickness. 
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FIG. 13: High resolution magnetic diffraction patterns from 
several SL samples measured with NG-1 reflectometer. The 
solid curves are fits of Eq. lO to the data points. The fitted 
values of the partial correlation coefficient p for each spectrum 
(as defined in subsection D) is shown in the figure. The small 
additional peaks visible in between the magnetic SL satellites 
must not be attributed to the correlations of opposite sign - 
they can as well result from different periodicity in a small 
portion of the SL (compare Appendix). 



B. Neutron diffraction in high magnetic field 



The evolution of the (^22) reflection for two, 
[(EuTe)4|(PbTe)i2]4oo and [(EuTe)5|(PbTe)i5]4oo, SLs in 
ffext parallel to the [110] axis is presented in Fig. 
With increasing i?ext; antiferromagnetically coupled 
spins flrst rotate towards directions perpendicular to the 
field and then gradually incHne towards the field (see 
schemes in Fig.EJ- Accordingly, the AFM diffrac- 
tion structure gradually fades away, while a new set 
of peaks emerges in the (111) FM diffraction region at 
Qx — ^- The almost total disappearance of the AFM 
component at 6 T is consistent with the behavior of bulk 
EuTe in external magnetic fieldsSi. With respect to the 
magnetic interlayer coupling, it is crucial to note, that 
the FM SL peaks at Qz = 1 are quite narrow (only 
sHghtly broadened beyond the instrumental Hnewidth) 
as compared to the significantly broader satellite peaks 
in the AFM region at Qz = 1/2 and Sext = 0- In ad- 
dition, there is no broad background in the FM region. 
The sharp satellite peaks in the FM state are due to the 
perfect long range spin coherence due to spin alignment 



by the high external magnetic field, which gives addi- 
tional clear evidence for the excellent structural quality 
of the samples. As a consequence, the broadening of the 
satellite peaks in the AFM region cannot be attributed 
to structural imperfections, but must be due to a limited 
long range spin coherency between the magnetic layers in 
the AFM state that may be induced only by spontaneous 
magnetic interlayer interactions. 



C. Field cooled samples 

All the above experiments have been carried out on 
samples which were cooled to temperatures below in 
zero external magnetic field. In this case, the magnetic 
field applied to a sample that is already in a correlated 
state does not destroy the existing interlayer correlation 
unless the field becomes strong enough to infiuence the 
AFM order within the individual layers. This takes place 
only for external fields much stronger than 1 T. In con- 
trast, cooling the samples from above to below in 
relatively weak magnetic fields (of the order of 100 - 200 
G) almost entirely prevents the formation of any inter- 
layer correlations. This is demonstrated in Fig. ll5[ where 
the magnetic diffraction patterns obtained for zero-field- 
cooled (ZFC) and the field-cooled (FC) superlattices are 
depicted. The magnetic field was applied parallel to the 
SL growth plane and oriented along the [110] crystallo- 
graphic axis. In the field-cooled cases (open symbols), 
the spectra only have the form of the structure factor of 
a single layer |-FblP, characteristic for uncorrelated SLs 
(see Fig. llir a)). whereas the spectra obtained after zero- 
field cooHng show the usual satellite peaks attributed to 
interlayer coupHng (as shown in Fig. [T]T &)). This change 
of the diffraction spectra just by the application of an ex- 
ternal field during cooling again demonstrates the purely 
magnetic origin of the multi-peak structure. Detailed ex- 
periments, in which the (5^^) peak profile was studied 
after cooling the sample in different external fields show 
that the loss of correlations exhibits a gradual depen- 
dence on the field strength. The multi-peak spectrum 
starts to evolve towards a broad maximum already for 
cooling at fields as low as 10 G, but the final uncorre- 
lated state is being reached only when cooling at higher 
fields, usually a few hundred gauss. However, the effect 
is fully reversible - subsequent warming up and cooling 
down the samples again in zero field restores the origi- 
nal correlated state. All the samples under investigation 
have shown this type of behavior. 

The possible explanation of the different behavior of 
the FC and ZFC samples will be given in the next sub- 
section. Here, however, it should be emphasized that 
such a behavior excludes the possibility that the satellite 
structure of the (^^1) peak results from the pinholes 
in the PbTe spacer and the formation of EuTe bridges 
across the spacer region, because in such case the result- 
ing interlayer coupHng would not depend so sensitively 
on the applied external fields. 



11 




0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 

coordinate (reciprocal lattice units) Qz coordinate (reciprocal lattice units) 



FIG. 14: Extended neutron diffraction scans for (4/12) (a) and (5/13) (6) SL samples, showing the suppression of the AFM 
scattering at the (^^|) reflection point (Qz = 0.5) with increasing external magnetic flelds, and the emerging of new satellite 
lines near the (111) reciprocal lattice point {Qz = 1.0) due to the induced FM spin alignment (chemical structure contributions 
are subtracted - the data represent purely magnetic scattering). The shaded bars (removed data points) indicate the regions 
where a very strong reflection from the BaF2 substrate occurs. The solid curves are the best flts of Eq. (O to the data. The 
corresponding alignments of magnetization in successive spin monolayers are shown by the arrows. 



D. Comparison of the experimentally determined 
interlayer correlations with theoretical predictions. 



Three mechanisms of the interlayer correlations, which 
can be relevant to AFM/nonmagnetic semiconductor 
layer structures, have been proposed in literature. First, 
the transfer of the magnetic order to the next magnetic 
layer via the spin-polarized carriers bound to the impuri- 
ties located in the spacer, was considered in Refs. 12^123 . 
This does not seem to apply to EuTe/PbTe SLs, since 
in PbTe the large dielectric constant and small carrier 
effective masses prevent the formation of shallow im- 
purity centers^. Another mechanism considered is the 
long-range dipolar interaction, which was investigated in 
R.ef ■ [2^ for FM metallic layer systems with domain struc- 
ture. It does not exist for perfect AFM layers, but one 
can argue that in real SLs the dipolar coupHng between 
local magnetic moments related to interface terraces and 
steps as invoked in Section llll (^1 can be effective. We ex- 



pect, however, this mechanism to be much weaker here 
than in the FM case, since the dipole-dipole interaction is 
proportional to the square of the average dipole moment. 
In addition, we expect that such mechanism should be 
much more effective for SLs with odd number of mag- 
netic monolayers than for those with m even, where the 
terraces do not lead to local dipole moments. No such 
preference was observed in the experimental data. 

Finally, in Ref. US a mechanism was presented, which 
attributes the interlayer coupHng to the sensitivity of the 
SL electronic energies to the magnetic order in consec- 
utive magnetic layers. The total energy of the valence 
electrons for two different magnetic SLs, one with the 
same and the other with opposite spin configurations in 
neighboring magnetic layers, was compared. The differ- 
ence between these two energies was considered as a mea- 
sure of the strength of the interlayer magnetic coupHng 
resulting from band structure effects. In Ref. it was 
shown that in both studied types of IV- VI semiconductor 
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FIG. 15: Neutron diffraction scans performed at (^^f) re- 
ciprocal lattice position for the zero-field-cooled (ZFC) and 
field-cooled (FC) (5/10), (5/12), and (4/12) SLs. The FC 
data has been rescaled to show the effect of enveloping the 
ZFC data. These diffraction data were collected on the NG-1 
refiectometer. The absence of the interlayer correlations in 
the samples after field cooling is evident (some traces are still 
visible in the 4/12 sample). For comparison the calculated 
single layer structure factor is plotted with a continuous line. 



structures, i.e., in FM EuS/PbS and AFM EuTe/PbTe, 
this mechanism can be eflFective. As shown in Fig.^J the 
calculated strength of the coupling decreases monotoni- 
cally with the increasing thickness n of the spacer and is 
practically independent of the thickness m of magnetic 
layers. In addition, it was found that for all FM and AFM 
SLs, regardless of their m and n values, the lower en- 
ergy corresponds to the antiparallel ahgnment of the two 
spins facing each other across the spacer layer. In other 
words, the mechanism leads to an antiferromagnetic cou- 
pling between the FM layers. In the AFM EuTe/PbTe 
structures, however, the energetically preferred spin con- 
figuration along the SL growth axis depends on the parity 
of the number of monolayers within the magnetic layer, 
i.e., the actual magnetic period is equal to the chemical 
period for even m but twice as large for odd m. 

These theoretical results have proven to explain the ex- 
perimental observations in the FM (001) EuS/PbS SLsli. 
For FM structures the magnetization and neutron refiec- 
tivity measurements in external magnetic fields, enable 
one to determine directly the strength of the coupHng 
and compare it with the model. The sign of the in- 
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FIG. 16: The interlayer exchange constant J, defined by the 
difference of the total electronic energy for the same and op- 
posite spin configurations in successive magnetic layers of the 
SL, as a function of the spacer thickness for FM EuS/PbS 
and AFM EuTe/PbTe SLs. Open circles represent the ex- 
perimentally obtained J for FM EuS/PbS. For the AFM SLs 
the direct experimental determination of the strength of the 
coupling is not possible. 



terlayer exchange coupling and the rate of its decrease 
with the PbS nonmagnetic spacer thickness are in very 
good agreement with the predictions of the model, as 
shown in Fig. El The fact that the experimental values 
of the exchange constants estimated from the saturation 
fields in real FM structures are about an order of mag- 
nitude smaller than the theoretical ones was attributed 
in Ref . ITll to the interfacial roughness and interdiffusion, 
which were shown to reduce significantly the strength of 
the interlayer coupHng, also in metaUic structures. 

For the AFM EuTe/PbTe structures the comparison 
of the theoretical predictions with the experimental data 
is much more complicated than for the FM structures 
- in this case not only the perfect tool to measure the 
strength of the interlayer coupling, i.e., the saturation 
magnetization, is not applicable, but, as shown below, 
the correlated spin configurations are much more sensi- 
tive to the morphology of the SL. 

In our superlattices evidence for the interlayer coupling 
between the AFM EuTe layers comes from the satellite 
structure of the neutron diffraction spectra. The sign of 
the coupling can be determined only by a detailed anal- 
ysis of the positions of the satellite peaks. Moreover, 
to describe the observed shapes of the AFM diffraction 
spectra, we have to invoke the idea of "partial correla- 
tions", described by an interlayer correlation parameter 
p {\p\ < 1), as presented in detail in the Appendix. An 
idealized fully correlated EuTe/PbTe SL would contain 
a single S'-domain, with the monolayer magnetization se- 
quence in any i^^ layer either repeated in the (i -I- 1)**^ 
layer (perfect correlations with p = +1) 

Titit Titit Titit Titit, 

or reversed in the (i-f l)**" layer (perfect correlations with 
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In both cases, the AFM neutron diffraction pattern 
should exhibit a series of very narrow peaks, with the 
width defined by the instrumental resolution only, as seen 
for the satellite peaks in the FM region at i/ext = 6 T 
in Fig.^^and Fig. lTTT bV However, in real samples, this 
perfect long range SL order may be disrupted by "phase 
lapses" (i.e., switches to the other S-domains types) oc- 
curring at random intervals. Such "partially correlated" 
chains can then be characterized by fractional p values, 
now expressing the probability P = (1 +p)/2 that any 
two adjacent EuTe layers have identical spin sequences. 
Applying diffraction theory to such a system, one obtains 
the following expression for the magnetic diffraction in- 
tensity (see the Appendix for the formula derivation): 



HQz) CC |i^BL(Q. 



l-2pcos{Q,D)+p^ 



(9) 



where D — m^EuTe + '^'5pbTe is the SL period; (JguTe 
and i5pbTe being the monolayer thicknesses of EuTe and 
PbTe, respectively. This expression is similar to that 
used for analyzing x-ray diffraction patterns from par- 
tially ordered layered structuresSS. The value of p deter- 
mines both the widths of the AFM satellite peaks as well 
as the height of the underlying "hump" (see Fig. ITTjl . 
By adjusting the p parameter for each sample the ob- 
served spectral shapes are reproduced remarkably well 
(soHd fines in Fig. [T3jl . The least-square fitted |p| val- 
ues are considerably lower than unity, even for relatively 
thin PbTe spacers. This indicates the presence of mech- 
anisms inhibiting the correlation formation in the AFM 
EuTe/PbTe SLs. Although one can speculate that sev- 
eral effects, such as thickness fiuctuations or interface 
roughness, may participate in the suppression, neutron 
diffraction measurements in moderate external in-plane 
fields -ffext < 1 T imply that magnetic anisotropy fields 
i/a in random directions parallel to the layers play a ma- 
jor role in this effect. Two ^-domains in EuTe layers 
facing each other across the nonmagnetic spacer may be- 
come correlated only if the interlayer coupfing energy is 
sufficiently high to overcome the anisotropy in at least 
one of them. 

This observation together with the model described 
above offers a simple qualitative explanation of the be- 
havior of the field cooled samples shown in Fig. El It is 
based on the fact that the interlayer coupHng becomes ef- 
fective only below the Neel temperature, when the AFM 
order in EuTe layers is already well established and the 
anisotropy fields are still weak. The interlayer coupHng 
energy resulting from band structure effects is propor- 
tional to the cosine of the angle 9 between the spins at 
the opposite borders of the nonmagnetic spacer. Thus, 
the torque responsible for the relative rotations of the 
spins in neighboring EuTe layers should be proportional 
to sin 0. During cooling of the sample in an external mag- 
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FIG. 17: Influence of the magnitude and the sign of the 
correlation parameter p on the magnetic diflraction intensity 
from 5/15 SL, as calculated from Eq. J^J. The change of the 
sign of p leads to the shift of the SL satellite peaks by half 
of the spectrum periodicity. With the decreasing value of \p\ 
the width of the peaks increases and so does the height of the 
underlying "hump". The single EuTe layer magnetic structure 
factor I-FblP, enveloping the whole spectrum, is presented 
with the dashed line. 



netic field, all the spins ahgn along the direction perpen- 
dicular to the magnetic field. Thus, the torque is equal 
to zero and the interlayer correlations cannot be formed. 
In contrast, in the ZFC samples the AFM ordered spins 
in different magnetic layers align randomly along differ- 
ent in-plane directions, hence sin6' and the torques are in 
general not equal to zero. Therefore, the rotation mech- 
anism can be effective, leading to the correlated layer 
structure in the case of ZFC samples. 

From Eq. ^ and Fig. El it follows that for a given 
combination of m and n values the change of the p sign 
corresponds to a half period shift ^AQsl in the AFM 
satellite positions. For Qz's in the the vicinity of (^^5) 
reciprocal lattice point the same is expected when the p 
is fixed, but either m or n is changed by ±1. This al- 
lows one to determine the relative spin configurations in 
successive layers in the SLs from the obtained sign of the 
interlayer coefficient. We note that this is possible under 
a strong assumption that the structures are morpholog- 
ically perfect, i.e., with the same, well defined m and n 
values throughout the entire (EuTe)m|(PbTe)„ SL com- 
posed of several hundreds of periods. From such an anal- 
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FIG. 18: Magnetic diffraction pattern from highly perfect 
EuTe/PbTe (5/15) x 10 SL (circulated solid line) showing the 
strong interlayer correlation. In this case the negative sign of 
p is in agreement with the theoretical model predictions. The 
experimental data for the less perfect SL with 300 repetitions 
(dashed line) is shown for comparison. 



ysis it turned out that the spectra for the superlattices 
with nominally even m and even n reveal the preference 
for the same monolayer spin sequence in successive EuTe 
layers (•• ti •• Ti ••)• F'or ^Ls with odd m and even n 
neutron diffraction spectra indicated reversed configura- 
tion (•• tit •• iti ••)■ Both these configurations are in 
agreement with the theoretical model prediction. In con- 
trast, for the samples with m and n both nominally odd, 
the neutron diffraction spectra seem to indicate that the 
(• • • tit • • • Tit • • •) configuration is preferred, contrary 
to the theoretical results. 

To shed light on this issue, an effort to detect the inter- 
layer exchange coupling in EuTe/PbTe SLs with a smaller 
number of SL periods, i.e., in SLs with better controlled 
m and n values, was undertaken. This task is not trivial 
as the intensity of neutron diffraction spectrum depends 
crucially on the number of spins involved. From the ad- 
ditional series of EuTe/PbTe SLs with only 10 periods, 
one sample indeed showed SL satellite peaks in positions 
corresponding to the sign of the coupling predicted by 
the model. In Fig. El we present the comparison of the 
spectra of two SLs, both with nominal m — 5 and n = 15, 
and with different number of SL periods TV. Clearly, the 
sample with only 10 periods shows the expected negative 
p value as compared to the opposite sign for the previ- 
ously measured sample with 300 periods (dashed line in 
Fig. CU. This result seems to suggest that in the long 
process of MBE growth of the SLs with a large number 
of repetitions (typical growth time of several hours) the 
preference occurs to form terraces with even number of 
monolayers. The reasons leading to such tendency re- 
main unclear, but one conceivable explanation may be 
that the number of monolayers which form a unit cell of 
the bulk material is somehow preferred during the long 
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FIG. 19: Dependence of the interlayer correlation parameter 
p on the PbTe spacer thickness in EuTe/PbTe superlattices as 
determined from the fitting of the neutron diffraction spectra. 



process of layer growth. In the [111] direction this cor- 
responds to two cation and two anion sheets. Such a 
mechanism would be valid for both, EuTe and PbTe, con- 
stituent materials. 

The theoretically predicted decrease of the strength of 
the coupling with the nonmagnetic spacer thickness is 
reflected in the decrease of the values of the correlation 
coefficient p, as shown in Fig.^l A quantitative compar- 
ison between the experimental and theoretical results is 
not possible in this respect. As far as the range of the in- 
teraction is concerned, the experimentally observed very 
long range of the interlayer interactions seems to exceed 
the range predicted by the model. The weak correla- 
tions still visible in samples with very thick spacers can 
be ascribed to the possible contribution from a residual 
dipole-dipole interaction. 



V. CONCLUSIONS 

We have performed an extensive study of the mag- 
netic properties of AFM EuTe epitaxial layers and of 
[(EuTe)m|(PbTe)„]Ar SLs grown by molecular beam epi- 
taxy along the [111] direction. The structural proper- 
ties of these samples were characterized by high resolu- 
tion x-ray diffraction. The magnetization and neutron 
diffraction experiments show that the magnetic proper- 
ties of the type II AFM EuTe layers depend sensitively 
on lattice distortions stemming from EuTe/PbTe lattice 
constant mismatch. Due to the resulting biaxial strain 
and finite EuTe layer thickness, the FM spin-sheets are 
all oriented parallel to the (111) growth plane, i.e., they 
form a single T-domain. 

For a large number of EuTe/PbTe SLs, a systematic 
study of the Neel temperature dependence on the number 
of EuTe and PbTe monolayers in the SL period was per- 
formed. The transition temperature to the AFM phase 
depends on the strain state of the magnetic layers as well 
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as on their finite thickness. The observed changes of the 
Neel temperature are described by the dependence of the 
exchange parameters on the lattice distortions and they 
follow essentially a mean field behavior. 

From neutron diffraction measurements in an applied 
magnetic field parallel to the (111) growth plane, detailed 
information on the S- type domain structure and on the 
in-plane anisotropy fields was obtained. The latter ones 
are considerably higher than in bulk EuTe. These fields 
play an important role in the formation of interlayer 
correlations. Neutron diffraction experiments in mod- 
erate magnetic fields and magnetization measurements 
showed that, irrespective of the number of monolayers in 
the EuTe layer, no net magnetic moment is present in 
the studied SLs. In this sense the EuTe/PbTe system 
constitutes a prototypical example of an antiferromag- 
netic/diamagnetic superlattice. 

The most interesting feature of these SLs is a pro- 
nounced interlayer spin correlation between successive 
EuTe layers revealed by magnetic neutron diffraction. 
The characteristic fingerprint of these correlations are SL 
satellite peaks in the vicinity of (HI) reciprocal lattice 
point. The correlations persist up to PbTe layer thick- 
nesses of about 60 A. Based on kinematical diffraction 
theory, the formula describing the diffracted beam inten- 
sity as a function of momentum transfer, Qz, has been 
derived for a general case of partially correlated SLs. A 
correlation parameter p (obtained by least-square-fitting 
to the neutron diffraction spectra) was found to follow 
a downward trend with increasing thickness of the non- 
magnetic spacer layer, thus, reflecting the weakening of 
the interlayer interactions with the distance between the 
magnetic layers. The signs of p - that govern the spin 
sequences in successive EuTe layers - were compared 
with the predictions of the theoretical model presented in 
Ref. lU In this model the interlayer correlations are me- 
diated by valence-band electrons and are inferred from 
the minimization of the total electronic energy of the 
EuTe/PbTe system on the spin arrangements in adja- 
cent magnetic layers. Essentially, the major features of 
this theoretical model, namely: 

(i) monotonic decay of the interlayer interactions with 
the distance between magnetic layers, 

(ii) independence of coupling strength on magnetic 
layer thickness, and 

(iii) opposite directions of the spins in the bounding 
monolayers of the two consecutive EuTe layers fac- 
ing each other across the PbTe spacer 

have been confirmed in our neutron and magnetization 
experiments, although in order to be able to check exper- 
imentally the last issue (especially for m — odd/n — odd 
SLs) samples with extreme structural perfection were 
necessary. 



APPENDIX: NEUTRON DIFFRACTION 
SPECTRUM FOR AFM SUPERLATTICES 

In this Appendix we calculate the profiles of diffraction 
spectra for a SL made up of alternating N antiferromag- 
netic layers and N nonmagnetic layers, each consisting 
of TO and n atomic monolayers, respectively. For simplic- 
ity, in the following it is assumed that in both SL con- 
stituent materials the spacing between the monolayers 
has the same value of d and that there are only two pos- 
sible directions of spin (magnetization) in each magnetic 
monolayer (it can be shown, however, that the results 
of this Appendix remain valid for systems with several 
in-plane easy axes and for layers consisting of many S- 
domains^^) . 

We consider three different situations: 

(a) perfectly correlated SLs - interlayer correlations 
lead to one of the two types of magnetic order, il- 
lustrated in the diagrams in Section IIV Dl in the 
entire SL; 

(b) uncorrelated SLs] and 

(c) partially correlated SLs - structures, in which there 
is a dominant tendency to form one type of corre- 
lations between the successive magnetic layers, but 
due to some disruptive mechanisms a minority of 
nearest-neighbor layers pairs is aligned in the op- 
posite way. 

In the standard kinematical theory approach the 
diffracted wave, resulting from magnetic scattering of un- 
polarized neutrons, is obtained by adding up all waves 
diffracted by individual magnetic atoms. This leads to 
the following equation: 

^diff PC f{Q) Kj exp{iQ ■ rj ) (A.l) 

3 

where Q is the scattering vector, f{Q) is the single-atom 
magnetic formfactor, fj is the position of the j^^ atom, 
and Kj is the magnetic scattering amplitude for a single 
atom, equal +k or —k for the "up" and "down" spin ori- 
entation, respectively. In the symmetric reflection geom- 
etry, most often used in diffraction studies of multilayers, 
the scattering vector is parallel to the superlattice axis: 
Q — (Q,0,Qz). The summation over individual atoms 
can be then replaced by a summation over the mono- 
layers. Since for all atoms located in the l*"^ monolayer 
Q ■ Tj — Qz ■ z — Qzld, the equation simplifies to: 

V'diff cx /(Q) J2 exp{iQzld) (A.2) 
I 

where Aii, the sum of magnetic scattering amplitudes of 
all atoms residing in the monolayer, is proportional 
to the monolayer magnetization. Taking advantage of 
the SL periodicity, one can separate this equation into a 
summation over all monolayers within a SL "elementary 
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cell" - a bilayer (BL) - and over the N SL repeats. Thus, 
for a bilayer, consisting of m magnetic monolayers (and 
n nonmagnetic ones, for which all the Mi — 0), one can 
define the magnetic structure factor Fbl as: 

m — 1 

Fbl{Qz) = f{Q) J2 cM^Q^l^d) (A.3) 
Eq. ljA.2|l can thus be written in the form: 

N-l 

V'diff OC 

where D = {m + n)d is the SL period. The spin configu- 
ration in the v^^ magnetic layer with respect to the first 
layer is described in Eq. ljA.4|l by the S,^ coefficient, which 
takes the value +1 for the same and —1 for the opposite 
magnetization sequences. 

The intensity I{Qz) of the diffracted radiation is given 
by: 

N-l N-l 

l^diffP ^J2Y1 ea?/3exp[zQ,(a-/3)]|FBL(g.)P 

(A.5) 

where the structure factor square |fBL(<3z)P can be writ- 
ten as: 



|i^BL(Q^)|' OC 



cos (mQ^rf/2) r , J 



(A.6) 



sin^(mQ^rf/2) 
cos2(Q^d/2) 



for m even 



The structure factor term \Fbl{Qz)\'^ has broad max- 
ima (with weak subsidiaries on both sides) centered at 



Qz = 5(27r/d), |(27r/(i), . . ., i.e., half-way in between 
the reciprocal lattice points corresponding to the basic 
atomic structure with periodicity d. 

Calculating the spectrum profiles I{Qz) for SLs 
with perfect interlayer correlations requires putting in 
Eq. IIA.511 £,1, appropriate for the given type of correla- 
tion. The task reduces then to summing geometric pro- 
gressions, which yields: for = I 



/(Qz)cx |Fbl(Q. 



and for £^ ^ (-1)''+^ 



HQz) (X |Fbl(Q.)P X 



^ sin\NQ,D/2) 
sin\Q,D/2) 



(A.7a) 



""'^fr?n7.f for iV even 



(A.7b) 

These functions consist of sharp maxima at regular 
intervals AQz — 2tt/D. The intensity of the narrow lines 
is "modulated" by the structure factor, what produces the 
characteristic groups of peaks. It should be noted that a 
change from one to another type of interlayer correlations 
causes the narrow line positions to shift by ^AQ^, which 
makes possible to detect such a transition. 

In order to analyze uncorrelated and partially corre- 
lated SLs the double sum in Eq. ljA.5|l should be re- 
arranged into sums over different kinds of layer pairs: 
namely, the sum of all same-layer terms (a = the 
sum of all terms with \a — (3\ = 1 (i.e., corresponding to 
adjacent magnetic layers), all terms with |a — /?| = 2 (i.e., 
corresponding to next-nearest layer pairs), an so on: 



N-l N-l 

E E ^"^^ exp[*Q,(a -(3)]^N+J2 ^^^P cos[QzD{a - (3)] 

Q=0 13=0 a=if3 
N-2 N-3 
N + 2 COsiQzD) + 2 C0s(2g^L>) CaCa+2 + ...+2 COs[{N - l)QzD]£_o£.N-l (A., 

I 



The £a£(3 product for any pair of layers, labeled a and 
/3, can be thought of as the correlation coefficient for this 
pair. The number of layer pairs that are kD apart is 
N — k; hence, the average correlation coefficient for all 
such pairs in the SL structure can be written as: 

^ N-k-l 



a=0 



Using Eqs. (|A.8ll and ljA.9|l one can write the Eq. ljA.5|l 



for diffracted intensity in a simple form: 

N-2 



I{Qz) (X iVl^^BLl' 



l + 2^Pfc(l- -)cos(Q,i?fc) 



k=l 



(A.IO) 

For a "perfectly random" superlattice the correlation 
coefficients for all layer pairs vanish on statistical aver- 
aging. Hence, for the uncorrelated system: 



I{Qz) (X N\Fbl{Q, 



(A.ll) 
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i.e., the diflraction spectrum reproduces the shape of the 
structure factor square - in agreement with the expected 
result for a random system with no coherence between 
the waves scattered by individual layers. 

The last situation to discuss is the partially correlated 
superlattice. If it is assumed that the only relevant in- 
teractions are between the nearest layers, and there are 
no long-range interactions, which introduce coupling be- 
tween the more distant layer pairs, then it is straightfor- 
ward to show that the correlation coefficient for second- 
nearest layers is p2 = Pi, for third-nearest layers is 
P3 — pf, etc. In the following we drop the subscript 
and denote pi by p. 

If the value of \p\ is significantly lower than 1, the cor- 
relation coefficients rapidly decrease, and only the first 
few terms in the sum in Eq. fOnli are relevant. If, in 
addition, the number of repeats N in the SL is large, one 
can use the approximation 1 — k/N = 1, and obtain: 



HQz) (X |Fbl(Q. 



l + 2^/cos(Q^i:'A:) 



A;=l 



(A.12) 



By applying the identity^fi: 



1 + 2^p''cos(fca;) = j 



l-p2 



fe=i 



— 2p cos(a;) + p'^ 



one obtains the final formula for the spectrum profile: 



HQz) « \FBLiQz 



1 -p2 



l-2pcosiQ,D)+p^ 



(A.13) 
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